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2 : Abstract 

I We consider momentum broadening and energy loss of high momentum partons in a hot non- 

■ Abehan plasma due to collisions. We solve the coupled system of Wong- Yang-Mills equations on 

^ . a lattice in real time, including binary hard elastic collisions among the partons. The collision 

^ I kernel is constructed such that the total collisional energy loss and momentum broadening are 



lattice spacing independent. We find that the transport coefficient q corresponding to transverse 
momentum broadening receives sizeable contributions from a power-law tail in the pj_-distribution 
of high-momentum partons. We establish the scaling of q and of dE/dx with density, temperature 
and energy in the weak-coupling regime. We also estimate the nuclear modification factor Raa 
due to elastic energy loss of a jet in a classical Yang-Mills field. 

PACS numbers: ll.15.Kc, 12.38.Mh, 24.10.Lx, 24.85. -hp, 25.75.Bh 
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I. INTRODUCTION 



The study of high transverse momentum jets produced in heavy-ion colhsions can provide 

tl 



information on the properties o 

iiflflflay.fly.fl 



region 



le hot QCD plasma produced in the central rapidity 
10 . 11 . Ql- After the discovery of jet quenching at the 



Relativistic Heavy Ion Collider (RHIC) 



13l . Il4i | a lot of progress has been made towards 



using jets as a quantitative tomographic probe of the QGP. Jet quenching refers to the 
suppression of high transverse momentum hadrons, such as vr*^ and rj mesons in central 
Au + Au collisions compared to expectations from measurements in p+p collisions. Whereas 
pions and ?7-mesons exhibit the same suppression at high p±, direct photons were found to be 



unsuppressed 



151 ]. This indicates that the observed suppression is related to the absorption 



(energy loss) of energetic partons in the medium. 

In this paper we study collisional energy loss and momentum broadening of massless 
high momentum partons traversing a non-Abelian plasma. Soft multi-particle interactions 
are treated by solving the coupled system of Wong- Yang-Mills equations in real time. In 
addition, particles can undergo hard elastic collisions. 

So far, estimates based on perturbative QCD (pQCD) of the strength of the coupling 
of jets to a plasma are sensitive to infrared cutoffs. Also, they are often restricted to 
systems that are (at least locally) in thermal equilibrium. The problem does not arise in 



the Wong- Yang-Mills simulation [16 



3, Q, Q, 



20 



21 



23 



24, 



25| . since the soft sector 



is described by classical chromo-fields. It is separated from the hard sector corresponding 
to hard elastic pQCD processes. The soft sector is non-perturbative but in an essentially 
classical way because the occupation number of field modes below the saturation momentum 

It is well known that a cutoff independent collisional energy loss can be obtained by 



30, 



3l|) we show 



resumming soft interactions [28|, |29|. In the present paper (see, also, Refs. 
by explicit implementation that this can also be achieved within the framework of a transport 
theory by treating the soft interactions via classical Yang-Mills fields defined on a lattice. 
We find that, in practice, this works even for physical values of the gauge coupling, ~ 2, 
so long as a weak-coupling (resp. continuum- limit) condition specified in Eq. (1251) below is 
satisfied. Within this framework, we are able to also consider the interesting problem of 
elastic energy loss of a jet propagating through a classical non-Abelian field, which might be 
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relevant for describing the early stages of a high-energy collision of large nuclei; see below. 

The main purpose of this paper is to extract lattice-spacing independent results for the 
transport coefficient q associated with broadening of the momentum distribution of gluon 



jets, as well as for collisional energy loss dE/dx. Previous publications 30|, l3l|] already 
presented a calculation of q within this approach, however lacking the detailed analysis 
shown here as well as a computation of dE/dx. Furthermore, in this paper we extract the 
entire p\ distribution of jets passing through a thermal plasma (not only its first moment q). 
We also address the scaling of q and dE/dx with the particle density, temperature, and jet 
energy. The scaling laws turn out to agree, qualitatively, with pQCD expectations although 
the overall magnitude of q and dE/dx is found to receive substantial corrections. 

We extract a value for q of 3.6 ± 0.3 GeV^fm"^ at T = 400 MeV in a thermal SU(3) 
background for a parton with energy E = 19.2 GeV. For the collisional energy loss we 
obtain dE/dx = 1.6 ± 0.4 GeV fm^^ 

This paper is organized as follows: We introduce the Boltzmann-Vlasov equations as 
well as the Wong equations for non-Abelian plasmas in Sec. [Tll and discuss the lattice 
implementation in Sec. Illli We outline how collisions are included into the Wong- Yang- 
Mills simulation in Sec. IIVI and describe how the separation between the soft and hard 
sector is done in Sec. |Vl After discussing the initialization of the simulation in Sec. IVIt we 
present results for collisional energy loss and for momentum broadening in Sec. IVIIi Finally, 
we close with conclusions in Sec. IVIIII 

II. BOLTZMANN-VLASOV EQUATION FOR NON-ABELIAN GAUGE THEO- 
RIES 

The classical transport theory for non-Abelian plasmas has been established by Heinz 



and Elze 



32 



33 



34 



35j. Here, we solve numerically the classical transport equation for 
hard gluons with adjoint SU(2) color charge q = q°'T°' , where the are the color generators, 
including hard binary collisions 

jf {d, + gq'^F^^^d; + gr"^A'^{x)q'^d,a) f = C. (1) 

/ = f{x,p,q) denotes the single-particle phase space distribution, F^^^ = d^A^ — d^A'^ — 

a 



gfo-bcj^b AC jg gauge field strength tensor, g the gauge coupling, A" the soft gauge field, 
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and C is the collision term to be defined below. It is coupled self-consistently to the Yang- 
Mills equation for the soft gluon fields, 



9 



dqqv" f{x,p,q) 



(2) 



with yf" = (l,p/p). For C = 



action near equilibrium 



36 



:hese equations reproduce the "hard thermal loop" effective 



37 



38|. However, the full classical transport theory (ITpi) also 



includes some higher n-point vertices of the dimensionally reduced effective action for static 
gluons 39(1 beyond the hard-loop approximation. The back-reaction of the long-wavelength 
fields on the hard particles ("bending" of their trajectories) is taken into account. This is 
essential for achieving cutoff independent results for the transport coefficient q and for the 
energy loss dE/dx of high momentum partons. 

When the phase-space density is parametrically small, / = which is the case for 

hard momenta, the collision term is given by 
1 



C 



»12 



/l/2|A^ 



12^1'2' 



, (3) 



P2 >^Pi "'P2 



with Jp = J (^2n)^2E- ■ '^^^ matrix element A4 includes all gg gg tree- level diagrams shown 
in Fig. [H and color factors as appropriate for the SU(2) gauge group. 
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O 
O 
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FIG. 1: Processes contributing to gg gg scattering at leading order. 



We employ the test particle method and replace the continuous distribution f{x,p, q) by 
a large number of test particles j4ol |: 



/(x,p,g) = ^Y.5\^~^m{'^^xf5^'\v-V^m^''"~'\q-q^{t)) 

^^test , 



(4) 



which leads to the Wong equations [16|] (also see j4l|, |42| ) 



Xi(t) = Vi(t) , 

Ut)=gqKt) (E"(t) + v,(t)xB'^(t)), 
qi{t) = -Wv':{t)[A^{t),q^{t)]. 



(5) 
(6) 
(7) 
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Here, Xj(t), Vj(t), and qi(t), are the position, velocity^ and color charge of the i^^ test particle. 
iVtest denotes the number of test particles per physical particle, A^{t) = y4^(t)r" and the 
commutator commutes color generators r". The last equation ([7]) describes the precession 
of the color charge qi{t) due to the color fields. 

Writing the current in terms of the individual test particles, the Yang-Mills equation for 
the soft gluon fields becomes 

= = ^ ^ q^^^Si^ - Mt)) . (8) 



test 



The theory without collisions as given by equations (jMZl) coupled to the lattice Yang-Mills 
equations ([8]) was first solved in [43] to study Chern-Simons number diffusion in non-Abelian 
gauge theories at finite temperature. It was applied later also to the problem of gauge-field 



instabilities in anisotropic SU(2) plasmas 2J,|25|]. Our numerical implementation is based on 



the improved formulation detailed in [25| where the non-Abelian currents, generated by the 
hard particle modes on the lattice sites, are "smeared" . This technique makes simulations 
in three dimensions on large lattices possible in practice. 

III. REAL-TIME LATTICE SIMULATION 

The time evolution of the Yang-Mills field is determined by the standard Hamiltonian 



method in = gauge [l7|, |43|, |4J] . The temporal gauge is particularly useful because it 



allows for a simple identification of the canonical momentum as the electric field 

= -A" . (9) 
In addition, time-like link variables U, defined below, become simple identity matrices. 



The lattice Hamiltonian in this gauge is given by 45] 



2 2 ^ Ntest,L 



including the particle contribution l/A^test,L \PLj\- The plaquette is defined by 

Un = U,my{t + x)Ul{t + y)Uli^ , (H) 



^ We consider only massless particles here so that |vi| = 1. 
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with the hnk variable 

Note that the index yU on [/ is merely an indicator of its direction and not a Lorentz index. 
The shifts x and y are one lattice spacing in length and directed into the x- or y- direction, 
respectively. ^ 

Eq. flTU]) is given in lattice units, which are chosen such that all lattice variables are 
dimensionless: 

with the lattice spacing a. Hi is hence related to the physical Hamiltonian by if = 
A/{g^a) Hl. To convert lattice variables to physical units we will fix the lattice length L in 
fm, which will then determine the physical scale for a. All other dimensionful quantities can 
then be determined from Eqs. (|T3l) . The Hamiltonian (ITOl) determines the energy density of 
the system and enters the equations of motion for the fields, e.g., 

jEi = {El,Hl}, (14) 

with the Poisson bracket {■,■}. Our lattice has periodic boundary conditions in all spatial 
directions. 



IV. COLLISIONS 



The col 



48 



49, 



50 



ision 



kernel 



51 



51, 



53 



54 



3l) is similar to that used in (parton) cascade simulations [46|, |47|, 



55l | . Here, it is restricted to hard binary collisions since soft multi- 



parton interactions are mediated by interactions with the collective Yang-Mills field. This 
way, we are able also to study collective phenomena and their contribution to isotropiza- 
tion and thermalization. In particular, we can in principle also study systems away from 
equilibrium (see 30|]) for which the scale corresponding to the Debye-mass squared in an 



isotropic system becomes negative 



56 



571, 



58 



59| . In this case it can obviously not damp 



the propagator to act as a cut off for the momentum exchange in the infrared. 



^ We set Ta = (Ta, the Pauli matrices, without the usual factor of 1/2, i.e., the commutation relation reads 
[T",r^] = 2(5°''. Another factor of 1/2 is absorbed into the ^-field, which has to be taken into account 
when calculating the physical fields E and B from it. 
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To complete our dual particle/field description, we need to specify the separation scale 
k* between the field and particle degrees of freedom. We will discuss this separation scale in 
detail below. For now it will serve as a lower bound for the exchanged momenta for binary 
elastic particle collisions. All softer momentum exchanges are mediated by the fields. 

The collision term is incorporated using the stochastic method introduced and applied 
601. l6l[| . We do not interpret the cross section in a geometrical way as done in 




52l | but determine scattering processes in a stochastic manner by sampling 



lision algorithm can be 



55l |. which will also be 



possible transitions in a volume element per time interval. This co 
extended to include inelastic processes gg ^ ggg as done in [54, 
incorporated in the future in our simulations. 

The collision rate in a spatial volume element A^x per unit phase space for a particle 
pair with momenta in the range (pi, pi + A^pi) and (p2, P2 + A.^P2) follows from Eq. ([3]) 

AiVcoU 1 A^p2 



Atj^A^xA^Pi 2Ei {2Tif2E2 



/1/2 



1 [ d^p'i d^p'2 



2 J {2Tif2E[ {2Txf2E^ 
Expressing the distribution functions as 

AAT, 



r|A<i2^i'2'|'(27r)^5W(pi +P2 -p'l -P^). (15) 



1,2, 



(16) 



and employing the usual definition of the cross section for massless particles 
1 f d^p[ d^p'2 



62| 



0"22 



-|-Mi2^i'2'r(27r)^(5W(pi +P2 -p'l -p'2) 



As J {27r)^2E[ (271)32^2 
one obtains the total collision probability in a volume element A^x and time interval At: 



(17) 



AN, 



coll 



At 



22 



frelO"22- 



(18) 



AN1AN2 A3x' 

■^rei = S/2E1E2 denotes the relative velocity, where s is the invariant mass of the particle pair. 
P22 is a number between and 1.^ Whether or not a collision occurs is sampled stochastically 
as follows: We compare P22 to a uniformly distributed random number between and 1. If 
the random number is less than P22, the collision does occur. Otherwise, there is no collision 



In practice one has to choose suitable A'^x and Ai to ensure that P22 < 1- 
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between the two particles in that time step. Since we represent each physical particle by 
iVtest test particles, we have to rescale the cross section as a ^ a/Ntest- This leads to 

To determine this probability, we require the total cross section a22- To leading order in 



it follows from the differential cross section obtained from the diagrams in Fig. [T] 63|, |6J, |65| : 

da Anal f tu su st \ , , 

with Nc the number of colors. The invariant Mandelstam variables are 



s = {pi+P2)\ t = {pi-p[Y, u = {pi-p'2f. (21) 

Using t = —q^i with q the momentum transfer, and the identity s + 1 + m = for massless 
particles, we can express the total cross section for processes with larger than k* as 

022= / -ndq\ 22 

The momentum transfer is then determined stochastically in the center-of-momentum frame 
of the two colliding particles from the probability distribution 

In Eq. f l22|) we have introduced the cutoff k* for point-like binary collisions. To avoid double- 
counting, this cutoff should be on the order of the hardest field mode that can be represented 
on the given lattice, k* ~ vr/a. 

V. SEPARATION SCALE 

The scattering processes in the regime of hard momentum exchange are described by 
elastic binary collisions, while soft momentum exchanges are mediated by the fields. A 
scattering in the soft regime corresponds to deflection of a particle in the field of the other (s). 

Physically, the separation scale k* should be sufficiently small so that the soft field modes 



below k* are highly occupied 4J] and hence can be described classically. On the other hand, 
k* should be sufficiently large to ensure that hard modes can be represented by particles and 
that collisions are described by ([3]), which is valid only for low occupation numbers since the 

8 



Bose enhancement factor (1 + /) is approximated by 1. In practice ~ 1, and we choose 
k* to be on the order of the temperature. At the same time, k* is related to the hardest 
available field mode, which on a cubic lattice is given by ^/Sn/a. Obviously, any matching 
between soft and hard regimes can only be done approximately, because the lattice on which 
the field modes are defined is cubic, while the momentum space cutoff of the hard collision 
integral is implicitly spherical. 
The "soft" scale is given by 



where A^c = 2 is the number of colors and n denotes the number density of hard gluons, 
summed over two helicities and A*"^ — 1 colors. Also, ph ~ 3T is the typical momentum of a 
hard particle from the medium. 

To allow for reliable numerical simulations one should have moL ^ vr and moCL ^ vr. 
The first condition ensures that the relevant soft modes actually fit on the lattice while the 
latter ensures that the lattice can resolve the wavelength I/itld to good precision. 

As we have argued above, we choose the inverse lattice spacing to be on the order of the 
temperature of the medium. Thus, with fl2^ the condition rriDa ^ tt roughly translates to 



In order to satisfy this relation, which is essentially the weak-coupling condition, at (7 ~ 1, we 
perform the numerical simulations below for an extremely hot and undersaturated medium: 
> n. This ensures that the simulations are carried out near the continuum limit. We 
verify below that transverse momentum broadening of a high-energy jet passing through 
a thermal medium is independent of T if the density and the ratio of jet momentum to 
temperature is fixed. One may therefore obtain a useful weak-coupling estimate of (p^) 
(resp. for the related transport coefficient q) by extrapolating our measurements down to 
temperatures relevant to present heavy-ion collisions. 

VI. INITIALIZATION 

We consider a heat-bath of Boltzmann distributed particles with a density of n = 
{5, 10, 20} fm"^ and an average particle momentum of 3T = {6, 12, 18, 24} GeV. For a given 




(24) 
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lattice (resp. k*) we take the initial energy density of the thermalized fields to be 



/d'^k 
j^;^ kfBose{k)&{k* - k) , (26) 



nTT^ 



where 

is a Bose distribution normalized to the assumed particle density n, and C, is the Riemann 
zeta function. 

The initial field amplitudes are sampled from a Gaussian distribution: {A'^{x)A^-{y)) = 
^6ij6°'^6{'x — y) . To thermalize the initial fields (approximately), we match their Fourier 
spectrum to the classical limit of the Bose distribution. Hence, the initial spectrum is 
gauge-fixed to Coulomb gauge and a filter is applied such that 

Ai ~ 1/k 

(in continuum notation). Setting Ei = initially"^. Gauss's law implies that the local charge 
density at time t = vanishes. We ensure that any particular initial condition satisfies 
exact local charge neutrality. The charge smearing algorithm for SU(2) explicitly exploits 
(covariant) current conservation and hence Gauss's law is satisfied exactly by construction 
at all times ^25]- 

The above procedure ensures that there is no large discontinuity of the energy density 
when going from the field to the particle regime. This way we are able to vary the separation 
scale k* about the temperature T by varying the lattice spacing. Fig.[2]shows the distribution 
of field modes and particles and the separation scale k* ~ T. 



VII. MOMENTUM DIFFUSION AND ENERGY LOSS OF HIGH MOMENTUM 
PARTONS 

Having initialized the background particles and fields, we can now add a few high- 
momentum test particles propagating along a given ("longitudinal") direction which rep- 
resent the jets. Their density should be sufficiently low so that they do not infiuence the 
thermal background significantly and so their mutual interaction is minimized. 

^ Equipartitioning of electric and magnetic fields is achieved very rapidly within a few time steps. 
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k/T 

FIG. 2: (Color online) Bose distribution and its low and high-momentum limits, used for the 
initial fields and particles, respectively. Physically, the separation k* should be on the order of the 
temperature T. The band between T/2 and 2T roughly indicates the region within which we vary 
k*. 



We always initialize "bunches" of test particles, which represent one physical hard mo- 
mentum parton ("jet"). A bunch corresponds to Attest particles in the same lattice cell. The 
physical color charge is independent of iVtest. If, in fact, the color charges of all test particles 
representing one jet add to zero, no coherent radiation is emitted (colorless jet). Such jets 
can only suffer coUisional energy loss^. In the particle-in-cell simulation radiative energy 



loss is not consistenly included (see e.g. 



Initializing a bunch of test particles with 



aligned color vectors, leading to a net current on the lattice, will hence not correspond to 
the correct physical bremsstrahlung process. We postpone the consistent implementation of 
radiative energy loss to future work. 

As detailed above, colorless bunches of test particles permit us to restrict to collisional 
energy loss and momentum broadening due to elastic collisions only. We first demonstrate 
that in our approach both 



^ Note that individual test particles from the bunch are of course colored and hence they collide not only 
with hard thermal particles but also with the modes of the thermal fields. 
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and the differential energy loss dE/dx are independent of the separation scale k* . Here and 
in what follows, p\_ denotes the momentum transverse to the initial jet momentum, q can 
be extracted from the squared transverse momentum of the test-particles accumulated up 
to a time t: 



Fig. [3] depicts the contributions to q due to soft and hard collisions, respectively, as well as 
the total. In these simulations the gluon density of the medium was taken to be n = 5 fm~^, 
the temperature T = 4 GeV, and the jet energy is 16 times the average thermal momentum 
(48 T). q is shown as a function of the separation scale k* ~ -\/37r/a. 

The curves for the total and the soft contributions were averaged over ~ 80 runs for 
each point. The curve corresponding to the hard sector was obtained by subtracting the 
result without hard collisions (soft sector only) from the total. The error bars indicate one 
standard deviation about the mean. 

We find that the total value is constant to a good approximation although the contribution 
due to soft scatterings changes considerably. Below k* ^ T, the soft sector contributes less 
than 10%, while it starts dominating around k* ^ 3T. It is evident, therefore, that transport 
coefficients obtained in the leading logarithmic (LL) approximation from the pure Boltzmann 
approach (without soft fields) are rather sensitive to the infrared cutoff k*, unless the energy 
^/s is extremely high. In LL approximation. 



where ~ s is the upper bound for the momentum transfer and C is a constant. 

Fig.Hlrepeats the same analysis for the collisional energy loss per unit path length, dE/dx. 
Again we find a constant total energy loss, and a similar dependence on k* of the partial 
contribution due to soft interactions as for q. We have also verified the /c*-independence 
for different temperatures, densities and jet energies. Thus, the above-mentioned matching 
of soft and hard processes provides estimates for q and dE/dx which are independent of 
the artificial separation scale k* (and of the lattice spacing a). It cures the infrared diver- 
gence of the perturbative hard-scattering cross section and does not rely on infrared cutoffs 
from equilibrium physics such as the Debye mass; thus, calculations are not restricted to 
equilibrium. On the other hand, the matching procedure might have to be modified for 




(30) 
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FIG. 4: (Color online) (In-)dependence of the energy loss dE/dx on the separation scale k* . T = 4 
GeV, g = 2,n = Sfni'^ 



other observables which are sensitive to very different scales. This should be analyzed in the 
future. 

Next, we turn to the density and temperature dependence of q and dE/dx. Figs. Oand 
inishow the linear rise of q and dE/dx with the density, which is expected from Eq. and 
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FIG. 5: (Color online) Linear dependence of q on the density n. T = A GeV, g = 2, E/T = 48, 
k* ~ 1.16 T. The line shows the best linear fit. 




FIG. 6: (Color online) Linear dependence of dE/dx on the density n. T = A GeV, = 2, E/T = 48, 
k* ~ 1.16 T. The line shows the best linear fit. 

the perturbative results (!30l) and (I3T1) below, respectively. We will use this linear dependence 
below to extrapolate to larger densities (e.g., ri^^^™'^^{T = 500 MeV) ~ 32fm~'^ for pure glue 
in SU(3)). 

Fig. [7] shows that q is approximately independent of T as long as the ratio of the jet 
energy to the temperature E/T as well as the density n are fixed. From (l30l) we expect 
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FIG. 7: (Color online) (In-)dependence of q on the temperature T. n = 5 fm ^, g = 2, E/T = 48, 
k* PS 1.16 T (squares) and k* « 1.73 7 (circles). 



at most a logarithmic dependence on T, because — s and (s) = GET. The simulation 
shows that this dependence is very weak. 

Fig. [8] shows dE/dx dropping approximately like ~ 1/T. This behavior is expected from 



the perturbative LL result for elastic energy loss (see, for example, ref. 67l |) 



-InU^'^^ , (31) 



dE _ / 16vra|Ag\ E f 

where s is the center-of-mass energy for a process involving scattering of the jet from a 
hard thermal excitation, ~ s is the upper limit for the momentum transfer, and C is a 
constant. Because (s) = GET this leads to dE/dx ~ 1/T. Additionally, T also appears in 
the logarithm, but this dependence turns out to be weak. 

We can now extrapolate to temperatures which are accessible in practice, for which direct 
computations can not be performed due to the numerical reasons explained above. For an 
ideal gas of thermal gluons at a temperature T the density n = 16T^({3)/tt'^ (for = 3). 
Using the linear dependence of g on n confirmed above, and its independence on T for 
fixed n and E/T, we find g ^ 7 ± 0.6 GeV^fm^^ at T = 500 MeV. This number has 
been rescaled to the color factors appropriate to SU(3) ^. Since E/T is fixed, this result 



^ We divide the results by the prefactors given in Eqs. ((30| and ((3T|l . respectively, which correspond to 
Nc — 2, and multiply by the prefactors appropriate for Nc = 3. 
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T[GeV] 

FIG. 8: (Color online) Dependence of dE/dx on the temperature T. n = 5 fm~'^, g = 2, k* ^ I.IGT 
(black squares) and k* ~ 1.73 T (red/grey squares), a . . . f are fit parameters - the plot shows the 
best fits. Extrapolation to temperatures T ~ 500 MeV leads to uncertainties of order 40 %. 



corresponds to a jet energy of = 48T = 24GeV. The quoted error arises from using 
different possible fits, including a logarithmic dependence on T or not, and from different 
choices for k*. For a temperature of T = 400 MeV and the corresponding thermal gluon 
density, we find q ^ 3.6±0.3 GeV^fm"^ (E = 19.2 GeV). We emphasize that our simulations 
do not account for quarks and anti-quarks which would provide a sizeable contribution to 
the thermal density. Nevertheless, such values for q are within the range extracted from 

y,y. 



RHIC data j68|. 

In Fig. [8] we present a possible extrapolation of dE/dx to temperatures around 500 MeV. 
We find dE/dx ^ 0.35 - 0.6 GeVfm"^ at T = 500 MeV, and for a jet energy oi E = 
48 T = 24 GeV. Adjusting the color factors as appropriate for SU(3) and extrapolating to 
the thermal density of gluons we find dE/dx ~ 2.5 ± 0.6 GeV fm~^. At T = 400 MeV, the 
result is dE/dx ^ 1.6 ± 0.4 GeV fm~i {E = 19.2 GeV). 

In Figs. [9] and [TO] we show how q and dE/dx depend on the energy E of the jet. The be- 
havior is logarithmic, in agreement with the perturbative expectation; compare to Eqs. (130]) 
and fl3ip . using ~ ET (the explicit factor of E in the numerator cancels since (s) = GET). 
A fit of the numerical result to ( l30i) . using ~ s, leads to C ~ 1.45, but is good only if 
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an additional prefactor of ~ 0.63 is allowed. This suggests that the perturbative result does 
not describe the numerical solution very well, which could perhaps be expected at g = 2. 
Repeating the analysis for dE/dx via Eq. (IHT]) . we find C ~ 7.7; again, a multiplicative 
factor needs to be included, this time it is ~ 0.14. Thus, for the jet energies considered here, 
there is a smaller "K factor" relative to pQCD at LL for dE/dx than for q; note that the 
former is sensitive also to longitudinal momentum exchanges while the latter is not. 

2 I , , , , 
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FIG. 9: (Color online) Jet-energy dependence of q for T = 4 GeV and n = 5 fm k* ^ 1.16 T. 
The line shows the fit to Eq. (j30p . with an overall multiplicative factor of 0.63. 



We have also determined the full ^^.'distribution of the high-momentum partons travers- 
ing the hot medium in order to assess the relative contributions from various processes to 
its first moment q. We find that over time the initial (5-function broadens to a Gaussian 
distribution with a power-law tail. This enhancement of transverse momentum broaden- 
ing reflects the well known result from QED and is in line with the findings for QCD in 
^72!. The enhancement has also been discussed within the higher twist formalism 



Refs. 

in Ref. 73|. What is perhaps less obvious is the relative magnitude of the Gaussian and 
power-law parts, which may be expected to be time dependent. However, for time scales 
typical of heavy-ion collisions we do not observe a large relative shift of these contributions. 

Figs. [TT] and [12] show the distribution of the high-momentum "jet" test particles after 
t ~ 2.8 fm and t ~ 5.2 fm, respectively, in a double-logarithmic plot versus p^/t. We scale p\ 
by the inverse time so that the basic features of the distribution are nearly time independent. 
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FIG. 10: (Color online) Jet-energy dependence of dE/dx for T = A GeV and n = 5 fm ^. k* 
1.16 r. The line shows the fit to Eq. ()3ip . with an overall multiplicative factor of 0.143. 
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FIG. 11: (Color online) ^^-distribution of the high- momentum {E/T = 48) "jet"-partons after 
t « 2.8 fm for T = 4 GeV and n = 20 fm^^ {q ^ 5.16 GeV^fm"^). 
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The low-px part follows a Gaussian distribution in p^. The power-law tail at large 
p_L behaves approximately as p^. This is expected for particles experiencing only few 
scatterings since in the high-energy limit the differential cross section da /dp\ ~ pj^, c.f. Eq. 
( l20l) . This is the probability distribution for the momentum transfer in a single hard collision. 
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FIG. 12: (Color online) ^^-distribution of the high-momentum {E/T = 48) "jet"-partons after 
t « 5.2 fm for T = 4GeV and n = 20 fm"^ {q ^ 5.16 GeV^fm"^). 



In both figures we also indicate the value of q to show that the power-law tail contributes 
significantly to this transport coefficient. For the densities, temperatures and jet energies 
considered here it is clearly not a very good approximation to determine the transport 
coefficient q from the Gaussian part of the distribution alone as this would underestimate q 
substantially: discarding the power-law tail from Figs. [TTIT21 gives gcauss ~ 0.6 GeV^fm"^. 

Note also that transverse momenta on the order of the temperature (T = 4GeV), such 
as the separation scale k* ~ 1.2T, correspond to p'^_i/t ^ 5.2GeV^fm~^ in Fig. [TT] and to 
pl/t ^ 3 GeV^fm"^ in Fig. [H Ab ove this value for p\/t the distribution is due almost 
entirely to hard collisions (we have checked that multiple soft collisions do not contribute 
much in that region). 

Finally, we also provide an estimate for the nuclear modification factor Raa of the jet 
spectrum due to elastic energy loss in a classical Yang-Mills field^. This is of relevance for 
collisions of heavy nuclei at high energies: the large number of gluons produced in the central 



rapidity region can be described as a classical field for a short time [l8|, [l9|, l20|, l2l|, |22|, |23 | 
until the field modes decohere and thermalize 25 1 . These classical fields produced in the 

n 

early stage of the collision also exhibit long-range correlations in rapidity which we 



^ Classical radiative energy loss has recently been considered in ref. ^74] but is not taken into account here. 
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presently neglect. 

We proceed as follows. From our simulations presented above, the elastic energy loss at 
a density n and for a separation scale k* (= -\/3f ) = 2T can be parameterized as 



dx iV2 - 1 6T V k*^ J 



(32) 



with K = 0.143 and C = 22.975. This form for dE/dx has been established numerically in 
the weak-coupling regime specified by Eq. fl25|) for A^^^ = 2; in what follows, we extrapolate 
it to Nc = ?) and to physical density and temperature. Note also that for such large k* (on 
the order of the so-called "saturation momentum" Qs) most of the energy density is due to 
the classical field. We evaluate this expression as a function of the jet energy for = 3, 
T = 400 MeV and the corresponding thermal density of gluons. This corresponds to an 
energy density of about 17 GeV/fm^, which is an app rop riate average over the first 1 fm/c 



of a central Au+Au collision at RHIC energy 
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FIG. 13: (Color online) Nuclear modification factor RAAip±) of jets due to elastic energy loss 
in a classical Yang-Mills field produced in the early stage of a relativistic heavy-ion collision at 
RHIC. The band indicates the uncertainty originating from the extrapolation of dE/dx to physical 
temperatures (compare to Fig. [8]). 



The initial transverse momentum distribution of jets at RHIC can be parameterized 
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approximately as 



' (33) 



with n A. Here, p_L denotes the momentum of a jet transverse to the colhding ion beams. 
In the central region [y ~ 0) it is equal to the jet energy. The final distribution due to 
interactions with the background is then given by 

dNf f .o , „/oN , , . . , . dNi 



1 dN, 



^ (l-ef. (34) 



(1 — e)2 d'^p'^dy 

Here, e denotes the fractional energy loss up to a time r, which we take to be 1 fm/c: 



e(P±) = -^(P±)- (35) 

p± dx 

Thus, the nuclear modification factor Raa at the parton level (neglecting hadronization) 
can be written as 

dNf/d^p^dy 

We find that €{p±) is on the order of 10% and that it decreases with increasing jet energy. 
However, due to the relatively steep initial spectrum of produced particles at RHIC, this 
can lead to ~ 30% - 50% suppression in the p_L-range between 5 GeV and 20 GeV; see 
Fig. [131 Clearly, the experimentally observed fiat Raa ~ 0.2 can not be accounted for fully 
by early-stage elastic energy loss in the classical field background. Nevertheless, our result 
shows that this contribution is significant and that it can not be neglected. 



VIII. SUMMARY AND CONCLUSIONS 



We have studied collisional energy loss as well as momentum broadening of high- 
momentum gluon jets in a hot and dense non-Abelian SU(2) plasma by solving the coupled 
system of Wong- Yang-Mills equations in real time on a lattice. We separate the soft from the 
hard momentum exchange interactions by introducing a separation scale k*. This separation 
scale is given by the inverse lattice spacing ~ 1/a, which determines the magnitude of the 
highest momentum field modes that can be represented on the lattice. We fix its physical 
value to be on the order of the temperature. Momentum exchanges below that scale are 
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mediated by the classical fields, those above the separation scale by direct elastic coUisions 
between the particles. The latter were implemented through the pQCD collision kernel and 
the stochastic method for determining scattering probabilities. 

We restricted to collisional energy loss by simulating effectively colorless jets (at the scale 
set by the lattice spacing). We were able to obtain lattice-spacing and hence separation- 
scale independent results for q and dE/dx in a static and weakly coupled plasma. The 
dependence on temperature and density, as well as on the jet energy was found to fol- 
low, qualitatively, expectations from pQCD. We then extrapolated our simulation results 
to thermal densities and to more realistic temperatures which could not be simulated di- 
rectly. For a thermal gluon plasma (no quarks and anti-quarks) at T = 400 MeV, and with 
color factors adjusted to the SU(3) gauge group, wc estimate q k, 3.6 ± 0.3 GeV^fm"^ and 
dE/dx 1.6 ± 0.4GeVfm~^, for a jet energy of 19.2 GeV. The errors are mainly due to 
the required extrapolation. At finite time (on the order of the transverse dimension of the 
collision zone in a heavy-ion collision), the p5_"distribution of the high-momentum partons 
is found to be well approximated by a Gaussian distribution at low p± and to exhibit a 
power-law tail ~ at high p±. The first moment of the distribution, i.e. the transport 
coefficient q, receives a large contribution from the power-law tail. 

We have also provided a first estimate of the (elastic) energy loss of a jet traversing 
a classical Yang-Mills field, which might emerge in the early stage of a collision of large 
nuclei at high energy. For a field energy density of about 15 GeV/fm^, the fractional energy 
loss over a time interval of r ~ 1 fm/c amounts to about 10% - 20%. Once convoluted 
with the steep ~ l/p± initial spectrum of produced hard particles, this results in a nuclear 
modification factor Raa — 0.5 - 0.8, indicating that energy loss in the early stage (before 
the onset of hydrodynamic behavior of the hot medium) may give a significant contribution 
to the observed Raa at RHIC. 
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